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Based on the dynamical mean field theory (DMFT) and angle resolved photoemission spectroscopy 
(ARPES), we have investigated the mechanism of high T c superconductivity in stoichiometric 
LiFeAs. The calculated spectrum is in excellent agreement with the observed ARPES measure- 
ment. The Fermi surface (FS) nesting, which is predicted in the conventional density functional 
theory method, is suppressed due to the orbital-dependent correlation effect with the DMFT method. 
We have shown that such marginal breakdown of the FS nesting is an essential condition to the 
spin-fluctuation mediated superconductivity, while the good FS nesting in NaFeAs induces a spin 
density wave ground state. Our results indicate that fully charge self-consistent description of the 
correlation effect is crucial in the description of the FS nesting-driven instabilities. 



PACS numbers: 74.70.Xa, 74.25.Jb, 75.10.Lp 

Iron pnictides have attracted much attention due to its 
high T c superconductivity (SC).[1] A prime candidate for 
the pairing glue is the fluctuating antiferromagnetic (AF) 
order. It has been argued that itinerant electrons form a 
spin density wave (SDW) via Fermi surface (FS) nesting, 
and that fluctuating SDW causes the SC transition in 
the vicinity of the SDW phase boundary. [2] In this sce- 
nario, a key ingredient to SC is the FS nesting property. 
However, the FS nesting property of iron pnictides has 
been controversial. Some systems are believed to possess 
good nesting properties, [3-5] while some are not. [6, 7] 
Also there have been many theories that emphasize the 
role of local Fe 3d electrons. [8-10] 

Resolving such controversy requires accurate determi- 
nation of FS topology and orbital characters. Theoretical 
simulations based on a first principles method can pro- 
vide such information. However, the conventional density 
functional theory (DFT) method often fails to describe 
the electronic structure due to the significant electron 
correlation effect in Fe-based superconductors. Calcu- 
lated bands have to be renormalized by approximately 
2~4 to fit the experimentally measured band width, and 
the predicted spin magnetic moment is about twice larger 
than the experimental value. On the other hand, the dy- 
namical mean field theory (DMFT) on top of the DFT 
showed consistent results with the measured bands, the 
anisotropy and the small magnetic moment. [11-16] 

Even though DMFT has been shown to work well, the 
comparison of calculated and measured band structures, 
especially the FS, is still important to check the valid- 
ity of the calculation. For an accurate comparison of 



theory and experiment, LiFeAs is the most suitable sys- 
tem at present. The most studied 122 (isostructural to 
BaFe2As2) systems do not have neutral cleavage planes, 
which affects surface sensitive techniques such as the an- 
gle resolved photoemission spectroscopy (ARPES). 

On the other hand, neutral cleavage surfaces of 
LiFeAs allow us an accurate experimental determination 
of the band structure, FS topology and quasi-particle 
dynamics. [6] In this respect, DMFT correction on the 
band structure of LiFeAs in comparison with experi- 
mental data should provide a unique opportunity to ob- 
tain accurate electronic structure information and resolve 
some of the issues on the origin of SC. 

For this reason, we performed both DMFT and 
ARPES studies on an intrinsic iron pnictide supercon- 
ductor LiFeAs. Our goal is to validate the accuracy 
of our DMFT method by comparing with ARPES re- 
sults, and to unravel the SC mechanism by analyzing the 
DMFT spectral function. It will be shown that the FS 
nesting is marginally suppressed in LiFeAs due to the se- 
lective correlation effect of each Fe 3d orbital, and that 
this gives rise to the electron pairing mediated by the 
spin- fluctuation (SF). 

Simulation on LiFeAs is based on the fully 
self-consistent DFT in combination with DMFT 
(DFT+DMFT) as implemented in WIEN2k. [17-20] The 
local self-energy due to the correlated Fe 3d orbital is ob- 
tained with the continuous time quantum Monte Carlo 
(CTQMC) impurity solver, where U=5.0 eV and J=0.8 
eV are used. [21] ARPES experiments were performed at 
HiSoR BL-9 and ALS BL-7 with similar conditions in 
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Ref.[22]. Single crystals used in the experiments were 
synthesized by Sn-fTux method. [23] 

In Figs. 1(a) and 1(b), the DFT band structures are 
shown by thin gray lines. Overlayed on top of the bands 
are the d XZjVZ (Fig. 1(a)) and d xy (Fig. 1(b)) orbital con- 
tributions, respectively, as indicated by the size of black 
circles. Two small hole bands near the Y and one small 
electron band near M are mostly from d xz ^ yz orbitals, 
while the large hole and electron bands are from d xy 
orbital. Using the self-energy E(u;) obtained from the 
DFT+DMFT calculation, we compute the momentum- 
resolved spectral function to inspect relative changes of 
those bands. The result is shown in the lower panel of 
Fig. 1(c). The DFT band structure in Fig. 1(c) is scaled 
by an average renormalization factor of 2.55 for compar- 
ison. The most noticeable aspect of the data near T is 
that two d XZj y Z related inner hole bands are selectively 
lowered, while a d xy related outer hole band is slightly 
shifted up, which will be discussed. 

The calculated spectrum can be compared with the 
ARPES data in the upper panel of Fig. 1(c) presented 
in the same scale. Three hole bands near T are shown 
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FIG. 1: (color online) DFT band structures (red lines) with 
(a) d xz , yz and (b) d xy contributions indicated by the size of 
black circles, (c) (upper panel) Second derivatives of ARPES 
data along the T-M direcition. (lower panel) DFT+DMFT 
spectral functions in red color scale map, and DFT band 
structures in green solid lines. DFT results are rescaled by an 
average renormlization factor of 2.55. 



clearly in both theoretical and experimental spectra. 
Near the M point, there is only one electron band in the 
ARPES, though the DFT+DMFT predicts two bands 
crossing Ep. The discrepancy can be ascribed to the 
matrix element effect. Indeed, one can observe in the 
experimental FSs in Figs. 2(a) and 2(b) that two bands 
cross the Ep near the M point where the orbital charac- 
ter varies with the Fermi wavevector. The vertical dashed 
lines in Fig. 1(c) denote positions of experimental Fermi 
wavevectors. Our theoretical predictions exhibit small 
deviation from experimental ones, especially for the outer 
electron pocket near M. This might be due to the usual 
Li deficiency in the samples. A better agreement is ob- 
tained if we lower the Fermi level by 0.01 eV as indi- 
cated by the horizontal dashed line in Fig. 1(c) (4% 
Li deficiency). [24] On the contrary, shift of the Fermi 
level alone does not give a good agreement between DFT 
(green sold lines) and experimental results. Especially, 
the hole bands deviate even further. 

The change of the FS size due to the correlation effect 
affects the nesting condition. Shown in Figs. 2(a) and 
2(b) are FSs at k z =0 and 7r, respectively. One can see 
that xz/yz-driven hole FSs shrink while the xy-driven 
hole FS expands compared to the DFT results. Similar 
results have been recently reported for various Fe-based 
superconductors. [12, 16] Our DFT+DMFT results also 
agree well with the experimental results on the right. 
Meanwhile, we do not see a sizable difference in the elec- 
tron FS size for the two methods, which results in su- 
pression of FS nesting between hole and electron pockets 
in the DFT+DMFT results. 

To check the nesting conditions, we show in Fig. 2(c) 
overlayed hole and electron FSs in the T-M-Z plane. 
While DFT result exhibits a good nesting between the 
second hole pocket and an electron pocket (left panel), 
nesting becomes poor for the DFT+DMFT result (right 
panel). Such difference in the nesting property can in- 
fluence the quantum transition behavior dramatically, as 
will be shown later. 

A natural question is on the mechanism of selective 
FS shift that increases (decreases) the xz/yz (xy) hole 
FS size and affects the nesting condition. We show 
that it is related to the charge transfer among Fe-3<i 
orbitals which is done during the fully self-consistent 
DFT+DMFT steps. We have checked the spectral func- 
tion calculated by the DMFT self energy without charge 
self consistency, and could not observe clear energy shifts 
as in the fully self-consistent case. Fig. 3(a) shows the 
variation of the DMFT charge density from the DFT re- 
sult. The electron density increases along Fe-Fe bonds 
(red color) whereas it decreases along Fe- As bonds (blue) . 
Since the cartesian x- and ?/-axes are chosen along Fe-Fe 
bonds in our calculation, the d xz j yz orbital occupancy is 
increased while the d xy orbital occupancy is decreased. 
This is consistent with the downshift of cL^/^-related 
bands and upshift of d^-related bands shown in Fig. 1. 
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Charge transfer from xy to xz/yz orbital can be under- 
stood to be caused by the difference in hybridization mag- 
nitude, t of each orbital. Under the common Coulomb 
interaction U, the renormalization factor, 1/Z oc U/t, 
varies from 2.1 for z 2 and x 2 — y 2 to 3.9 for xy. xz/yz 
has an intermediate value of 2.9. Therefore, the xy or- 
bital experiences a smaller hybridization (or larger local- 
ization) compared to xz/yz. It means a larger Coulomb 
energy cost, which favors less occupation in xy. Since ti g 
orbitals mainly contribute to near Ep states, xy electrons 
will be transfered to xz/yz orbitals. 

The charge redistribution due to the electron correla- 
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FIG. 2: (color online) DMFT (solid lines) and DFT (dashed) 
Fermi surfaces at (a) /c z =0, (b) k z = tt are compared with the 
experimental FSs. There are three hole pockets (hi, h2, h3) 
near Y and two electron pockets (el, e2) near M. The small- 
est hole pocket (hi) has strong k z dependence and does not 
appeat at k z = tt. (c) k z dependence of FSs along the Y-M-Z 
cut. To check the nesting condition, electron pockets denoted 
by dot-dashed lines are shifted by (7r,7r) and plotted at Y on 
top of hole pockets (solid lines). DFT and DFT+DMFT re- 
sults are shown in left and right panels, respectively. 



tion also substantially contributes to bonding properties. 
The decrease in the electron density along the Fe-As bond 
direction makes the Fe-As bond weaker, which induces a 
longer Fe-As bond or higher arsenic height compared to 
the DFT estimation. In Fig. 3(b), we show the calcu- 
lated total energy as a function of As height, h(As), by 
DFT and DFT+DMFT methods. Compared to the ex- 
perimental value of h(As) e£C p=1.5 A, [25] the DFT method 
predicts a position 0.1 A lower, which indicates overesti- 
mation of the binding in the DFT method. Meanwhile, 
the DFT+DMFT estimation is in good agreement with 
h(As) ea jp. 

As shown above, our DFT+DMFT results are in 
good agreement with the experimental measurements 
and shows that the nesting is marginally established. 
This suggests that LiFeAs is located near the SDW 
boundary and that the SF can be the pairing glue for the 
SC. In order to check such possibility, we have performed 
quantitative analysis of the nesting-driven instabilities by 
utilizing the calculated spectral function A(k, e). In Fig. 
4(a), we compare the bare susceptibilities Xo(q)- The 
nesting enhancement at q = M is much suppressed by 
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FIG. 3: (color online) (a) Isosurface plot of electron density 
difference between DFT and DFT+DMFT results. The red 
(bule) color means increase (decrease) of 0.003 e/A 3 upon the 
DMFT calculation. The cartesian x or y axis is chosen along 
the Fe-Fe bonds, so the increase is related mainly to d xz or d yz 
orbital of Fe, while the decrease to d xy . (b) Calculated DFT 
and DFT+DMFT total energies against As height relative to 
the Fe plane. 
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FIG. 4: (color online) (a) Bare susceptibilities Xo(q) calcu- 
lated by using the DMFT and DFT methods, normalized by 
Xo(q = 0). (b) Intra-orbital scattering contribution to Xo(q)- 
(c) and (d) SDW and SC instability plots versus temperature 
for LiFeAs (Lilll) and NaFeAs (Nalll), respectively. 



including the correlation effect. It is mainly due to the 
size mismatch between hole and electron FSs as shown 
in Fig. 2(c). Since the A(k, e) is mostly from xz/yz and 
xy orbitals, we focus on their contributions to Xo(q)- I n 
addition, only the intra-orbital scattering is considered 
because it gives the dominant contribution to the SC. [26] 
In Fig. 4 (b), the xz/yz intra-orbital contribution is sig- 
nificantly larger (about twice) than that of xy. It is be- 
cause there are more xz/yz orbital derived states near 
the Fermi level. The shape of Xo(q) m Fig. 4(a) is also 
well reproduced by the scattering between xz/yz domi- 
nant states. Therefore, the FS nesting can be understood 
to be effectively between two xz, yz derived bands (hole 
and electron bands near T and M, respectively). Note 
that the smallest hole pocket at T has a negligible effect 
on the nesting property. 

Suppressed nesting predicted by the DMFT spectrum 
implies that the SDW phase is also suppressed and the 
SF driven SC ordering becomes more probable. To have 
a SDW phase, the following condition is satisfied for a 
given nesting vector Q. 

KWo(Q,T) (1) 
A similar condition follows for the SF-driven SC: 



(2) 



where Vsf(<1) = Vg DW xo(q) / [1 -Vsdw Xo(<l)] denotes 
the singlet interaction channel, and Xsf(<1) is the as- 
sociated SF susceptibility. Using the calculated DMFT 
xz/yz spectra, we compute [1 — VsdwXo] and xsf at the 



nesting vector Q = M as a function of temperature in 
order to check the leading order between SC and SDW. 
Fig. 4(c) shows the critical Tsdw with the condition 
in Eq. (1). Also Tsc is estimated from the function of 
Xsf(Q) using the condition of Eq. (2). Note that Xsf(q) 
is a smooth monotonic function and does not need to be 
positive as shown in Fig. 4(c). Since Vsf(q) shows a sin- 
gular behavior near the SDW ordering for given Q and 
^SDW, the condition in Eq. (2) can be always realized 
when Xsf(c\) is positive. So Tsc can be defined where 
Vsf(Q) becomes positive as decreasing temperature. In 
LiFeAs, Tsc is always estimated to be higher than Tsdw 
as shown in Fig. 4 (c). 

For comparison, we show in Fig. 4(d) the SC and 
SDW instabilities of NaFeAs which has a SDW ground 
state. One can see that Tsdw is higher than Tsc m 
NaFeAs, in agreement with experiments. Such difference 
between NaFeAs and LiFeAs comes from the fact that the 
FS nesting is stronger in NaFeAs (in the DMFT calcula- 
tion), which enhances the SDW instability and increases 
the Tsdw- On the other hand, SF is suppressed due to 
the stable SDW phase and thus Tsc is decreased. As a 
result, our DFT+DMFT method correctly predicts both 
SC and SDW ground states of LiFeAs and NaFeAs, re- 
spectively. Note that the GGA spectrum does not give a 
SC transition because both compounds show good nest- 
ing features, resulting in a stable SDW ground state. 

In summary we have shown that the inclusion of elec- 
tron correlation effect is essential to correctly describe 
the orbital occupation, FS sizes, and As height, all 
which are important factors for unraveling the SC pairing 
mechanism. [27] Our fully self-consistent DFT+DMFT 
calculation has accurately described the marginal FS 
nesting of LiFeAs observed in ARPES experiment. Us- 
ing the DFT+DMFT spectrum, we also have successfully 
reproduced the SF-mediated SC in LiFeAs and SDW 
ground state in NaFeAs. 
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AN EFFECTIVE TWO-BAND MODEL HAMILTONIAN 



One of two bands is a hole band centred at (0,0) and the other is an electron band at 
(7r,7r). Those are driven by Fe d xz / yz orbital, which gives rise to intra-orbital scattering The 
effective Hamiltonian is given below. 

H = H + H A + H M (1) 

H = J2 tc(k)<L<*a + E 0(P)4a/ P « (2) 

k,Q; p,Q! 

^ = E < + &, < 4 » 

(5) 

£ c (k) and £/(p) represent dispersions of hole band c and electron band /, where k is relative 
to (0,0) and p is relative to (7r,7r). 

The SC pairing interaction Ha is characterized by the coupling constant V sc . 

H ^ = \Y^ ^*%?'a'( k >P) ( C L C -k/3/-P/3'/p«' + /L/-k/? C -P/3' C pa') 
k,p 

^W(k, P) = V^^U(ta% lal (7) 

(8) 

The SDW interaction H M is characterized by the coupling constant V SDW . 

h m = -\ E V X'(P'P; kk') (fl, a c P ^,f kW + fi p , a c- P ecl kp ,f_ kW ) (9) 

p'-p=k / -k 

Kw'a'KP P> kk ) = Vp;kk'^ ' ( 10 ) 



(6) 



SDW INSTABILITY 

We ignore the full momenta dependence of V SDW in Eq. (S10). So V^p.^, is assumed as 
a constant Vsdw- Electron-hole scattering SDW susceptibility within RPA can be written 
like below. 

( \ Vsdw / in x 

XSD "' (q) = 1 - WXb(q) <U) 



(<*) = i E / * / *' A /(p + * e )^(P' 60 /(e |_^ (£/) (f 2) 



The spectral functions A/(p, e), A c (p, e) are computed from our DFT or DFT+DMFT 
results. From obtained Green's function 6r#(p, e), where z,j represents band index, the 
whole hole spectrum is computed by ( — 1/tt) $^<i ImG^(p, e). Likewise the whole electron 
spectrum is (—1/tt) J2i>i c I m ^n(P> e )- Here z c is a chosen band number which separates hole 
and electron spectra. Finally obtained hole and electron spectra are projected into xz/yz 
orbital of Fe. Using this procedure, we show computed A/(p, e) and -A c (p, e) in Fig. SI 
along r~M direction near the Fermi level. One can see clearly the hole and electron bands 
near V and M, respectively. Also Af and A c behave like quasiparticle, so they are applicable 
to the effective two-band model. We note that one inner hole band among two xz/yz-drived 
bands in Fig. Sla has minor effect in drawing our main conclusion. 

At T = Tsdw the following condition is satisfied for a given nesting vector Q. 

l = V SD wXo(Q,T = T SDW ) (13) 




SPIN-FLUCTUATION MEDIATED SUPERCONDUCTIVITY 



The singlet interaction channel is given by 

V SF (q, i«) = X(q, <0) - *W = ^Wo(q,i») (M) 

i - K S wXo(q,Mi) 

Assuming s± pairing symmetry, the SF-driven superconducting instability is derived using 
associated hole (electron) Green's function G c (/) • It will be applicable to interacting spectral 
functions by DFT+DMFT. 

1 = \i Yl Yl G ^ p ' iuj )V S F(q, ityG c (p + q, iu + iO) 

iuj,ifl p,q 

xG/(-p, -icc;)^ j p(q, zfi)G c (-p + q, -iu + zfi) (15) 
The right hand side can be transformed in a following form involving A c ^y 

I de I de ' I dei I de[[V SF (q,in)] 2 A / (p,e)A c (p + q,e / )A/(-p,e 1 )A c (-p + q,e / 1 ) 

^ P iQ p,q ^ ^ ^ ^ 

1 1 1 1 

^ e — iu) e f — iu) — iQ €i + iu e\ + iu — iQ 
Also the last term is reexpressed by using Fermi-Dirac function f(e). 



1^ 1 1 



1 1 1 



B ^ e — iu e' — iu — iQ ei + iu e\ + iu — iQ ifl + e — e' iVt + ei — e\ 
/(e) + /( £l ) f(e>) + /( Cl ) /(e) + /ft) + /(e') + /ft) 



(17) 



e + ei -itt + e' + ei -iQ + e + e[ -2iQ + e' + e[ _ 

By using Vsf(<1, ity ~ Vsjr(q, zfi = 0), thus Vgi^q, zfi) — >■ VsF(q), we finally obtain a 
following form 

1 = ]t±YJ2 j de J de> J dei f de[{Vs F (c[)} 2 A f (p,e)A c (p + q,e')A f (-p,e 1 )A c (-p + q,e' 1 



) 



p q 



1 1 17(e) + /(e') + /(e!) /(e) + /&) | f(e') + f(e[) (ig) 



e — e' ei — [ e + ei e 7 + ei e + e[ e f + e[ 

We simplify the above expression as below. 

u iElMq)] 2 X5F(q) (19) 

q 

Computation of Eq. (S18) involves sixth-order integration including BZ summation, 
which is very challenging. For the fourth order integration over the energy, we use normal 
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Monte-Carlo sampling method. Important parameters that we need to care for good ac- 
curacy are (i) the number of p or q points, (ii) the number of MC sampling points for a 
given energy integration range, i.e. energy window, \e — /i| < E w . The BZ integration over 
p or q is done with small number of points as possible. So we choose the minimal grid size 
along z-direction that can give reasonable results, which is two. Also it will be seen that 
the summation over q is not necessary to judge the SC transition above SDW temperature. 
The E w is in principle to be infinite. But we choose such E w large enough to include most 
spectra of specific hole and electron bands. So this parameter is dependent on a problem. 

THE EXISTENCE OF T c ABOVE T S dw 

Noting that V&p(q) is singular at the SDW nesting vector Q, we express Eq. S19 as 
below . 

= ^ [Vsc(Q)] 2 Xsc(Q) + ^ £ [^c(q)] 2 Xscid) (20) 

Here N is the number of q points in the full Brillouin zone, which is supposed to be finite 
number. 

We suppose following two conditions, (i) Xsc(<l) > only at Q, (ii) dxsc J^ ,T " } < 0. The 
second condition is naturally satisfied in order that physically meaningfull SC transtion be 
observed at low temperature. Also in Eq. S18, the FD function ^+^^ satisfies the second 
condition. 

In Eq. S20, the second term at the right hand side is a negative quantity. At the Tsdw, 
if Xsc(Q) > 0, singular Vsc(Q) is large enough to satisfy the following condition. 

[Vsc(Q)] 2 Xsc(Q) > Nf3n 4 - (N - 1) £ [^c(q)f Xsc^) (21) 

Thus, from the condition (ii), there exists a SC transition which satisfies the Eq. S19 
above T^l^- This is also true for the case that Xsc(<l) involves positive quantities partially 
at some q points. 
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